Counting statistics for mesoscopic conductors with internal degrees of freedom 
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We consider the transport of electrons passing through a mesoscopic device possessing internal 
dynamical quantum degrees of freedom. The mutual interaction between the system and the con- 
duction electrons contributes to the current fluctuations, which we describe in terms of full counting 
statistics. We identify conditions where this discriminates coherent from incoherent internal dynam- 
ics, and also identify and illustrate conditions under which the device acts to dynamically bunch 
transmitted or reflected electrons, thereby generating super-Poissonian noise. 
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The probabilistic nature of quantum transport consti- 
tutes a fundamental source of current fluctuations, whose 
study provides detailed information about the carrier dy- 
namics [l|. A comprehensive characterization of these 
fluctuations is provided by the framework of full count- 
ing statistics (FCS), which delivers a unifying perspective 
on the conductance G = I/V (where / is the current and 
V is the bias voltage), the shot noise P — 2eIF (where 
F is the Fano factor, with F = 1 signifying the Pois- 
sonian statistics of uncorrelated carriers), and higher- 
order stationary current correlations 2] . In the regime of 
single-particle electronic transport close to equilibrium, 
the Pauli principle fundamentally constrains the current 
fluctuations, which results in sub-Poissonian noise with 
Fano factor F < 1. Repulsive interactions between the 
charge carriers generally result in a further reduction 
of the noise 0. On the other hand, super-Poissonian 
statistics can be achieved in strongly correlated systems, 
where interactions enable the formation of a charge den- 
sity wave 0| ■ The noise can also be enhanced via external 
means, e^g., by the presence of classical multi-state fluc- 
tuators [5(. 

Many quantum conductors of present interest, such 
as molecules with different conformations and nanoelec- 
tromechanical systems, possess internal quantum degrees 
of freedom which can be associated to multiple internal 
states of a mesoscopic device [see Fig. (Ha)]. These de- 
grees of freedom are not directly associated to the charge 
carriers, but cannot be considered external because of 
the backaction of the passing carriers. The purpose of 
this work is to demonstrate that these internal dynam- 
ics can be effectively probed via their influence on the 
counting statistics. We find that even though the internal 
degrees of freedom do not introduce strong correlations 
in the conventional sense |f|, they can result in super- 
Poissonian statistics with F ^> 1. The transport then 
takes the form of dynamically created trains of transmit- 
ted or reflected electrons. We identify conditions of ideal 
dynamical bunching and show that these can be realized 
in a model system, a coherent which-path interferome- 
ter 0, H| consisting of an Aharonov-Bohm ring which is 
electrostatically coupled to an excess electron in a dou- 



FIG. 1: (a) Sketch of a mesoscopic device with internal de- 
grees of freedom, attached to leads contacted by electronic 
reservoirs, (b) Realization of such a system in terms of a 
coherent which-path interferometer, consisting of an excess 
electron in a double quantum dot (states | t), | 1), tunnel 
splitting A), which interacts via a state-dependent potential 
4 with charge carriers passing through an Aharonov-Bohm 
ring. 



ble quantum dot [see Fig. []Tb), and results in Fig. [2], 
We also identify a regime where the transmitted bunches 
consist of exactly two electrons, with F = 2. 

These phenomena are identified by extension of pre- 
vious works on counting statistics in passive systems, in 
particular, counting statistics of conductors coupled to 
quantum detectors (electrometers) 0, [l(| , and the wave- 
packet approach [ll|, which provides a convenient dy- 
namical perspective on FCS. The resulting expression for 
FCS takes the form of a generalized propagation of the 
density matrix for the internal degrees of freedom. A key 
conceptual ingredient appearing in our expressions is the 
extent of coherence maintained throughout the evolution 
of the internal system dynamics, which manifests itself 
in electronic interference contributions associated to dif- 
ferent dynamical histories of the internal system state. 

Model and framework. — We consider electronic trans- 
port through a mesoscopic system with internal degrees 
of freedom, as schematically depicted in Fig. HJa). The 
system possesses N internal states \n), n = 1, 2, 3, . . . , N, 
which affect the transport of electrons passing through 
the system. In terms of creation and annihilation opera- 
tors i/j^ , ip of electrons in momentum (k) and position (r) 
representation, a reference Hamiltonian representative of 
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this problem takes the form 



wave-packet variant of the scattering approach, intro- 



nm 

+ J2 I dv^{v)i,{v)V n {v)\n){n\. (1) 



Here, the first term describes the kinetic energy of elec- 
trons with dispersion relation e(k), which is independent 
of the internal system state. The second term is the 
Hamiltonian of the isolated internal system degrees of 
freedom. The last term describes the interaction between 
these degrees of freedom and the electrons, which here is 
formulated in terms of a potential energy V n whose state- 
dependent part is assumed to be confined to the scatter- 
ing region (the potentials V n may also contain a state- 
independent part outside this region, which can model, 
e.g., some features of the leads). 

Note that we do not assume that the isolated system 
Hamiltonian commutes with the interaction part. There- 
fore, even when the internal degrees of freedom are ini- 
tially prepared in an eigenstate \a) of H sys (correspond- 
ing to energy E a , where a = 1, 2, 3, . . . , N), the interac- 
tion with scattered electrons will induce dynamics into 
the internal system. Moreover, since the system is not 
driven externally, total energy conservation implies that 
any change of internal system energy will be compen- 
sated by a change of energy of the transmitted or reflected 
electrons, which amounts to inelastic (but still coherent) 
scattering. We are interested in the consequence of these 
phenomena for the average current and its temporal fluc- 
tuations. 

A general framework for the statistical description 
of the electronic transport is provided by full counting 
statistics, which delivers information about the proba- 
bility P(Q; T) that a number Q of charge carriers are 
transmitted during a time interval T. This information 
is encoded into the moment generating function 



;(A) = J2 P(Q;T)ex P (XQ) = J^iQ* 



Q=0 



k=0 



(2) 



where (Q k ) denotes the fcth moment of Q. Of particular 
interest are the associated cumulants of order m, 



((Q m )) = 



d T ' 



<9A" 



ln X (A) 



(3) 
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which for large times increase linearly in T, ((Q m )) ~ 
(T/T)q m . Here 1/r is the attempt frequency of trans- 
mission for incoming electrons (r = h/eV for bias volt- 
age V"); the coefficient qi — g = G/Go then delivers the 
dimensionless conductance (in units of the conductance 
quantum Go — e 2 /h), while (72/91 = F delivers the Fano 
factor 

General results. — In order to obtain a general expres- 
sion for the moment generating function, we adopt the 



duced in absence of internal dynamics in Ref. 11 . The 
formalism takes a particularly simple form if we assume 
that at most one electron at a time is interacting with the 
internal system. This corresponds either to weak inter- 
actions, or to a small attempt frequency, where the latter 
can be realized either by applying a small bias voltage, or 
by judiciously injecting sparse charge carriers via an elec- 



tronic pumping device [12|. Exploiting that under these 



conditions wave packets of consecutive electrons are non- 
overlapping, the initial state of J\f electrons incident from 
the left lead can be formulated as the product of wave 
packets [l3| 



*imtiai({r„};t) = 



E^»> ( 4 ) 



where p a are the initial probability amplitudes of the 
internal system. Here, t n ~ nr denotes the arrival times 
of the electrons, which are spread out according to the 
attempt frequency, and Af ~ T/t. The final state of the 
reflected or transmitted electrons then takes the form 



*fi„ai({r„};t) =J2 



71=1 



where 



(5) 



(6) 



consists of a reflected and a transmitted wave packet, 
confined to the left and right lead, respectively; the in- 
dices indicate the final and initial internal state. These 
wave packets correspond to asymptotic transmission and 
reflection probabilities 

R*>* = J dr|rv CT (r;t>0)| 2 , (7) 
T ctV = J dr|^ ff (r;i»0)| 2 . (8) 

Following the general strategy of counting statistics de- 
veloped for systems without internal degrees of freedom, 
we formulate the moment generating function in terms of 
a generalized wave function, x(A) = (*final(A)|*fi n al(A)>, 
where 'I'finai(A) is of the same form as Eq. ([5]) but with 
^,(° ut ) replaced by 

^°fM,A) =r a , G {v ] t)+e x /H ala {v-t). (9) 

Here, the counting field A serves for book-keeping of the 
transmission events. After integration over the electron 
coordinates, the generating function takes the compact 
form of a matrix product, 



X (A) = X T M^ P , 



(10) 
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where 

Mv't'^s = J drrl, a (r;t)r ala (r;t) 

+ e A J ' dvtl, d {v;t)t al(J (v-t) (11) 

is the superoperator propagating the density matrix p 
of the internal system state, generalized to include the 
counting field A. Note that M. can be interpreted as 
a matrix where each row and column is specified by two 
indices, with one index each arising from the propagation 
of the bra and the ket of the internal system state (p then 
corresponds to a vector). For the pure state specified in 
Eq. the initial density matrix is p a & — p a p%] with 
modification of this initial condition, the formalism also 
applies if the initial internal state is mixed. The vector 
X with components = &a'a> embodies the trace of 

the final density matrix. In this dynamical description, 
the matrix elements M a i >C! = M. a i a i Ml7 are associated to 
the evolution of populations (the diagonal elements of 
p), while the other matrix elements describe corrections 
due to coherence (interference of electronic wave packets 
generated by different histories of the internal dynamics). 

In practice, the internal degree of freedom will suffer 
from decoherence due to coupling to the environment. 
This can be described by a modified generalized propa- 
gator M. = A4Z, where the diagonal matrix Z embodies 
suppression of off-diagonal elements of the intermediate 
density matrices. If the internal system state is contin- 
uously monitored, the internal dynamics become inco- 
herent, Zo-'er'.a-a- = &o> a> ^Ga^aa' , and the time evolution 
reduces to a classical stochastic problem described by a 
rate equation, which only involves transition probabili- 
ties between populations. The generating function then 
simplifies to 

X (A) = Y T M"P, (12) 

where 

M a . a = R a , a + e x T a , a (13) 

is the classical propagator of the population probabili- 
ties with initial conditions P CT = \po-\ 2 , while the vector 
Y T = (1,1,1,...) with unit components embodies the 
summation over all final states of the system [1J] . 

Of particular interest is the steady state behaviour of 
the system, corresponding to T > r. To analyze this 
limit we diagonalize the matrices M. (for coherent inter- 
nal dynamics), M (for incoherent internal dynamics), or 
M. (for cases in between these two limiting scenarios), 
and denote the largest eigenvalue by -D max (A). For large 
times, the cumulant generating function then becomes 

ln X (A) ~ -InZWA). (14) 

T 

Therefore, the problem of full counting statistics is re- 
duced to finding the eigenvalues of the matrices Ai, M, or 



M, respectively. Expanding D max (X) = Y^h=o^k^ k /k\ 
as a Taylor series, the dimensionless conductance is given 
by g = . d\ , and the Fano factor is given by F — (d 2 — 
d\)/d 1 M- 

Dynamical electron bunching. — In the conventional 
case of phase-coherent electronic transport through pas- 
sive devices, shot noise is suppressed due to the Pauli ex- 
clusion principle, resulting in a Fano factor F < 1 which 
indicates sub-Poissonian noise statistics P, 0- As we 
now show, even though we consider electrons that never 
directly interact with each other, the internal system de- 
grees of freedom can induce dynamical correlations which 
result in super- Poissonian noise with F> 1, correspond- 
ing to dynamical bunching of electrons. This possibility 
is most transparently revealed in the case of incoherent 
dynamics of an internal two-state system. In anticipa- 
tion of the concrete example discussed below, we denote 
the two internal states as a = +, — . The matrix M then 
takes the form 

M _(R ++ + e x T ++ + e*T + _ \ 

M -\R-+ + e x T |_ + e x T ) ' (15) 

and the largest eigenvalue can be calculated explicitly. 
Consider now the case of small switching probabilities 
S+- = R+- + T+_ < 1, S-+ = R-+ + T_+ < 1. On 
average, the two states are then populated with proba- 
bilities P+ = S+-/(S+- + S-+) and P_ = S-+/(S+- + 

S |_). In terms of these probabilities, the dimensionless 

conductance takes the simple form g ps P + T ++ + P_T__ 
of a time-averaged transmission probability. Assuming 
T ++ 7^ T , the Fano factor 

F „ 1 2P + p_(r ++ -r__) 2 

~ S+-+S-+ P+T++ + P-T— ' V ' 

on the other hand, can become arbitrarily large. The 
noise is then dominated by the rare probabilistic switch- 
ing between states of different transmission probabilities. 
In the limiting case that one of the two transmission prob- 
abilities is much larger than the other (say, T ++ 3> T__), 
the device is partially transmitting during some time in- 
tervals, but effectively shut down during others, which 
results in well spaced-out trains of transmitted electrons. 
These observations also extend to multi-state dynamics; 
noise is generally enhanced whenever the transmission 
probabilities differ between the states of the device. 

Application to a mesoscopic device. — We now illus- 
trate how conditions of dynamically enhanced current 
fluctuations can be realized in a specific mesoscopic de- 
vice, the coherent which-path interferometer shown in 
Fig. QJb) @, ■ The internal degree of freedom arises 
from an excess electron in a double quantum dot, with 
states | f)i I I) f° r occupation of the upper or lower 
dot, respectively. Depending on its location, this elec- 
tron blocks the path of charge carriers moving through 
the upper or lower arm of an Aharonov-Bohm ring, where 
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FIG. 2: Fano factor for the coherent which-path interferome- 
ter sketched in Fig. _0[b), where (a) shows results as function 
of the magnetic flux $ penetrating the Aharonov-Bohm ring, 
for fixed tunnel splitting A = 0.01 Ef, and (b) shows results 
as function of A, for fixed $ = h/2e. Solid (dashed) lines: co- 
herent (incoherent) internal dynamics. Insets in (a): dimen- 
sionless conductance (top), ratio of incoherent and coherent 
Fano factors (bottom). The result are obtained by propaga- 
tion of minimal-uncertainty wave packets with energy spread 
AE — 0.007 Ef- The interaction Vf t \. is modeled by an ef- 
fectively impenetrable barrier of height 20Ef, which blocks 
one of the two arms (of length L, with h 2 /2mL 2 — 0.061-ZJf) 
depending on the occupation of the double dot. 



coherent transport can be tuned by varying the mag- 
netic flux 3>. We assume that the isolated system eigen- 
states are the symmetric and antisymmetric orbital states 
l^) = (I t) i I l))/v / 2 of the excess electron, and that 
they arc separated by a tunnel splitting energy A (cor- 
responding to £± = =pA/2). 

To characterize the transport through this device, we 
calculate the matrix elements (|11[) by numerical prop- 
agation of initially Gaussian wave packets. The con- 
ductance and Fano factor are then obtained from the 
eigenvalues of the matrices Ai(X) (for coherent inter- 
nal dynamics) and M(X) (for incoherent internal dy- 
namics), as described above. Figure __[a) shows the re- 
sulting Fano factor as a function of $ for a set of sys- 
tem parameters. We observe that the current fluctu- 
ations can be significantly enhanced above the thresh- 
old F = 1 for super-Poissonian noise. For $ — > 0, the 

switching probabilities S-\ , S |- vanish because of the 

parity symmetry of the system [_J, and the Fano fac- 
tor increases over all bounds, as predicted by Eq. (|16| . 
In this regime, the Fano factor discriminates coherent 
from incoherent dynamics (bottom inset; as shown in the 
top inset, the conductance is unaffected by this distinc- 
tion). In contrast, when half a flux quantum penetrates 
the ring ($ = h/2e), parity symmetry ensures that each 
transmission event is strictly correlated with a transi- 
tion of the internal system degree of freedom, while no 
such transitions occur when an electron is reflected (i.e., 
= = T ++ = T__ = 0) 8]. Even in this 



case, dynamical electron bunching can occur: The maxi- 
mal Fano factor for incoherent internal dynamics now is 

F = 2, which is realized when R ,T_ + < i? ++ ,T + _ 

(or vice versa), so that P + ^ P_. The transport then 
consists of pairs of transmitted electrons, which coincide 
with the rare transitions into the state |— ) (first trans- 
mitted electron), followed by its immediate depopulation 
(second transmitted electron) . Figure __[b) shows that 
such conditions can be met by tuning the tunnel splitting 
A, and is insensitive to the degree of coherence in the in- 
ternal dynamics. Guided by these results, we expect that 
super-Poissonian noise statistics should be common when 
mesoscopic systems possess internal degrees of freedom. 

Concluding remarks. — In summary, we have investi- 
gated the phase coherent transport through mesoscopic 
devices with internal degrees of freedom, and formulated 
the full counting statistics in terms of a generalized prop- 
agator of the density matrix for these freedoms [Equa- 
tions (fT0|) -([T4 l) ]. In general, the counting statistics de- 
pend on the degree of coherence of the internal system 
dynamics. The dynamical switching of the internal sys- 
tem degrees of freedom can enhance the current fluctua- 
tions, resulting in Fano factors that significantly exceed 
the threshold F — 1 for super-Poissonian noise. 
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